Five loop Konishi from AdS/CFT 



Zoltan Bajnok"0, Arpad Hegedus^§, Rcjmuald A. Janik'^§ 
and Tomasz Lukowski^ 



° Theoretical Physics Research Group, Hungarian Academy of Sciences, 
1117 Budapest, Pazmany s. 1/A, Hungary 
''Research Institute for Particle and Nuclear Physics, Hungarian Academy of Sciences, 
H-1525 Budapest 114, POB 49, Hungary 
° Institute of Physics, Jagiellonian University, 
ul. Reymonta 4, 30-059 Krakow, Poland 



Abstract 

We derive the perturbative five loop anomalous dimension of the Konishi 
operator in N = 4 SYM theory from the integrable string sigma model by 
evaluating finite size effects using Liischer formulas adapted to multimagnon 
states at weak coupling. In addition, we derive the five loop wrapping con- 
tribution for the L = 2 single impurity state in the (3 deformed theory, which 
may be within reach of a direct perturbative computation. The Konishi 
expression exhibits two new features - a modification of Asymptotic Bethe 
Ansatz quantization and sensitiveness to an infinite set of coefficients of the 
BES/BHL dressing phase. The result satisfies nontrivial self-consistency con- 
ditions - simple transcendentality structure and cancellation of fi-teim poles. 
It may be a testing ground for the proposed AdS/CFT TBA systems. 
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1 Introduction 



The AdS/CFT correspondence [Ij states the equivalence of type IIB super- 
strings on AdS^ X 5*^ with the maximally supersymmetric Af = 4 gauge theory 
in four dimensions. Local operators in gauge theory correspond to states in 
superstring theory. In the planar limit, A^^c ^ oo, multitrace operators fac- 
torize essentially into independent single trace factors which correspond, on 
the string side, to noninteracting strings. In this limit, one identifies the 
spectrum of a type IIB superstring in AdS^, x 5*^ with the anomalous di- 
mensions of single trace operators in the maximally supersymmetric A/" = 4 
four-dimensional SU{Nc) gauge theory. 

This identification should hold for any value of the 't Hooft coupling con- 
stant A = Qym^c- In the strong coupling limit, and for operators with suffi- 
ciently large charges, the corresponding strings may become classical/semi- 
classical and then the problem becomes tractable by conventional methods. 
In contrast, once the coupling is no longer large and/or we consider generic 
short operators the answer requires us to consider the worldsheet QFT of 
the superstring on the quantum level. In particular, if we would like to make 
contact with conventional perturbative gauge theory calculations and explic- 
itly test the correspondence or follow some specific operator all the way from 
weak to strong coupling we would have to quantize exactly the worldsheet 
QFT of the superstring which is a highly nonlinear theory. 

Fortunately, the worldsheet QFT of the superstring in the AdS^ x 
background is integrable, which was first shown on the classical level in [2], 
and assuming that integrability also holds on the quantum level (for which 
there are now many indications), one may use the theory of two-dimensional 
integrable quantum field theories to eventually quantize the theory in a quite 
explicit way for any value of A. The AdS/CFT correspondence thus allows 
us to use tools and methods which exist only for two-dimensional quantum 
field theories to study the four- dimensional Af = 4 gauge theory. 

The theory of two-dimensional integrable quantum field theories is well- 
developed by now. The general strategy to determine their finite volume 
spectrum goes as follows: First one considers the model in infinite volume. 
The Hilbert space of asymptotic states is built up from noninteracting mul- 
tiparticle states which transform covariantly under the global symmetry al- 
gebra. Time evolution is formulated in terms of the scattering matrix that 
connects the initial and final multiparticle states [3]. Once integrability is 
assumed, there is no particle creation, and moreover any multiparticle scat- 
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tering process is shown to factorize into pairwise two particle scatterings. So 
the whole scattering information is contained in the 2 — 2 particle S-matrix. 
The requirements of unitarity, crossing symmetry, invariance under the global 
symmetry and the Yang-Baxter equation usually determines the scattering 
matrix uniquely up to CDD type ambiguities [4j. These are fixed through 
an analysis of the singularities of the scattering matrix all of which must 
have a physical origin. Thus the appearing poles have to correspond either 
to bound-states or to Coleman-Thun diagrams (anomalous thresholds). In 
this bootstrap solution we consider the bound-states and the original parti- 
cles on equal footing and determine the scattering matrix of the bound-states 
from the scatterings on their individual constituents [5j. Then the singularity 
structure of the bound-state scattering matrix is analyzed and new bound- 
states are searched for. The so-called bootstrap program is completed if all 
singularities of all the scattering matrices are explained and then the theory 
is completely solved in infinite volume. 

In the context of applications to the AdS/CFT correspondence, the infi- 
nite volume solution is not enough since one wants to describe closed strings 
and therefore consider the integrable quantum field theory on a cylinder of a 
given circumference related to a charge of the given state. In contrast 

to conventional relativistic theories, all (integer 0) sizes of the cylinder are 
in fact relevant for the complete spectrum of the theory. 

The finite volume solution of the model can be achieved by systematically 
taking into account the finite size effects due to the scatterings of particles. 
The leading finite size effect of a multiparticle state comes from the quanti- 
zation of momenta. It is described by Bethe-Yang equations and takes into 
account the scattering matrix in determining the allowed momenta |7j. It 
incorporates all polynomial corrections in the inverse of the volume. In ad- 
dition, there are exponentially small (Liischer) corrections as well and their 
leading contributions come from the polarization of the sea of virtual parti- 
cles |8|. For small volumes, these effects become dominant and one needs to 
perform a resummation of the virtual corrections, which sometimes can be 
carried out in the form of nonlinear integral equations. 

Indeed, the exact description of the finite volume ground state energy 
can be obtained from the Thermodynamic Bethe Ansatz (TBA) [9j. This 
is based on the fact that the contribution of the ground state energy domi- 

^See also [6] where it was shown that the meromorphicity of Y-system leads to the 
quantization of the temperature of the mirror model. 
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nates the Euclidean partition function for large Euclidean time, and that the 
same partition function can be calculated by exchanging the role of space 
and time. One is left with the determination of the partition function at 
finite temperature, but in the large volume limit where finite size effects are 
under control. This method provides typically coupled integral equations for 
pseudo energies which determine the ground state energy exactly. In some 
circumstances a careful analytical continuation (in the volume say) can pro- 
vide integral equations for excited states as well [10]. The procedure is a 
numerical one and the resulting equations are only conjectural which have 
to be further tested (but in all known cases they have passed all checks). 

The analogous bootstrap solution of the AdS^ x worldsheet QFT of 
the light-cone quantized Green-Schwartz superstring is currently almost com- 
pleted. Historically, the developments which led to it concentrated on the 
Bethe Ansatz part mostly on the gauge theory side |TTl [12], |13l HH [15] as 
well as for classical string solutions in AdS^ x [16], [17]. Later this was 
reformulated in the (spin-chain) S-matrix language in [181 [IH] 5 culminating 
in a proposal for the all loop Asymptotic Bethe Ansatz [TH] the derivation 
of the (spin-chain) S-matrix from global symmetry properties [SO]- In fact, 
this derivation could be transformed in a verbatim way into a starting point 
for the bootstrap solution of the worldsheet QFT. This was necessary since 
the worldsheet QFT perspective was crucial in order to tackle a new 'topo- 
logical' class of Feynman graphs - 'wrapping diagrams' - which went beyond 
the Bethe Ansatz in the spin chain guisqj. 

Let us therefore recall the basic steps of the bootstrap solution of the 
light-cone quantized Green-Schwartz string in AdS^ x 5*^. It is classically 
integrable for any value of the light-cone momentum, P+ identified with the 
f/(l)i^ charge J, which serves as the volume of the two-dimensional theory. A 
notable new feature of this theory is that it is not relativistic invariant. In the 
decompactification limit J — oo the massive excitation transform under the 
global symmetry algebra: the centrally extended sMc(2|2)^. This symmetry 
|20j . together with unitarity and crossing symmetry [22] completely fixes the 
scattering matriji0 including the dressing factor [211 [25] modulo ODD ambi- 
guities [2S]- The analysis of the pole structure revealed an infinite tower of 
bound-states [27| whose scattering matrices have been calculated as well [2^ . 
Double poles corresponding to anomalous thresholds were identified in [29] 

^Note however the attempt to use the Hubbard model formaUsm for this purpose |21j . 
^Taking into account some important subtleties |23| . 
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and the related Coleman-Thun diagrams were found. As the physical region 
of the AdS S-matrix is not known the complete analysis of its singularity 
structure is rigorously not completed yet. 

In order to find the spectrum one has to pass to finite volume and consider 
the theory on a cylinder. The corrections to the Asymptotic Bethe Ansatz 
can be identified with wrapping corrections [SU] , and the leading Liischer cor- 
rections for single [21] and multiparticle states [52] have been derived for this 
(nonrelativistic) theory. In this case the space-time interchange necessary for 
formulating TBA leads to quite a different theory [33j with a distinct set of 
bound-states. By now TBA equations are also developed by various groups 
[Ml [3S], under various analyticity assumptions, in particular about the 
analytically continued dressing phase. There are conjectures also for excited 
states TBA equations [35] but they have not been tested yet beyond the 
leading perturbative order of wrapping corrections. 

The formalism of the leading multiparticle Liischer corrections have been 
quantitatively tested in the case of the 4-loop anomalous dimension of the 
Konishi operator where the leading order wrapping diagrams have been cal- 
culated perturbatively [371 ES]- The corresponding leading order Liischer 
calculation found an excellent agreement [39]. Subsequently wrapping inter- 
actions computed from Liischer corrections were found to be crucial for the 
agreement of some structural properties of twist two operators [39] with LO 
and NLO BFKL expectations [40j. 

The aim of the present paper is to elaborate the Liischer correction to 
next to leading order and by this to calculate the anomalous dimension of the 
Konishi operator at five loops. Besides this explicit knowledge, which can 
serve as a testing ground for excited state TBA equations, we further test 
several issues of the formalism as well as a sizeable part of the BES/BHL 
dressing phase since the calculation requires the knowledge of the analyt- 
ically continued dressing phase in the Liischer kinematics. It relies on the 
conjectured finite size energy correction of multiparticle states, which at sub- 
leading order contains corrections originating from the modification of ABA. 
Although there is no perturbative gauge theory computation so far, the in- 
ternal consistency of our calculation provides enough confirmation to believe 
in its correctness. In addition we consider the five loop subleading wrapping 
correction for a single impurity operator in the f5 deformed theory, which 
may be within reach of a direct perturbative verification. 

The paper is organized as follows: In section 2 we summarize the main 
features of the 5-loop Konishi computation emphasizing the new phenom- 
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ena which appear w.r.t. the previous 4-loop case, and the possible internal 
consistency checks of the computation. In section 3 we rederive in a simple 
example the formalism of multiparticle Liischer correction. We focus on a 
theory where the S-matrix does not depend on the difference of the rapidities 
and derive TBA equation for the ground state. Excited state energy levels 
are obtained by analytical continuation and, by analyzing their large volume 
asymptotics, multiparticle Liischer corrections are extracted. As they con- 
tain the analytically continued scattering matrix, in Section 4 we determine 
the analytical continuation of the dressing phase by two different methods. 
Section 5 is devoted to the main computation of the anomalous dimension 
of the Konishi operator. It starts with the calculation of the ABA. Then 
we turn to the computation of the Liischer correction. It has two sources: 
one comes from the modification of the ABA, the other corresponds to the 
virtual particles circulating around the cylinder. Both terms include an in- 
tegration over the momentum of the mirror particles and a summation over 
their spectrum. Integration is carried out by residues, where, in comparison 
to the 4-loop case, we have to take into account an infinite tower of poles 
coming from the polygamma function part of the integrand. The resulting 
expression is composed of rational and polygamma functions which have to 
be summed over the bound-states. The technique developed for the summa- 
tion is explained in Section 6. Finally we give our conclusions in Section 7. 
The paper is followed by two Appendices. In the first we calculate the 5-loop 
anomalous dimension of a single impurity operator, which acquires nontrivial 
wrapping corrections in the (3 deformed theory. In the second Appendix we 
explain how to sum up terms containing polygamma functions. 

2 Main features of the 5-loop Konishi com- 
putation 

The Konishi operator tr $f is the simplest operator not protected by super- 
symmetry which, thanks to its short size has proved to be a testing ground 
for the AdS/CFT correspondence. In most computations it is more conve- 
nient to use a different representative of the same supermultiplet which lies 
in the si (2) sector 

tY{DZDZ)-tT{ZD'^Z) (1) 



6 




Figure 1: F-term of the Liischer correction corresponds to virtual particles 
propagating around the worldsheet and scattering with the Konishi two- 
particle state, a and h denote internal states of the virtual particle while q is 
its 'mirror' momentum. 

Its anomalous dimension following from ABA is given bj{§ 

Eaba = 4 + 12/ -48^^ + 336/ - (2820 + 288C(3))^^ (2) 

+ (26508 + 4320C(3) + 2880C(5))^^° + . . . 

However already at four loops there appears a contribution of wrapping 
graphs. The four loop wrapping contribution was computed directly in per- 
turbative gauge theory [37J using supergraph techniques and reconfirmed 
together with the nonwrapping part using component Feynman graphs in 
|38j . In [21] the same result was computed from the string sigma model in 
AdS^ X and came from a Liischer type F-term graph (see fig. 1) where a 
'virtual' particle was circulating in a loop. 

In the Liischer corrections the virtual particle is strictly speaking on-shell, 
however its kinematics are from the space-time interchanged theory (so-called 
mirror theory). One can estimate its contribution to the magnitude of the 

^See section Em for a quick derivation. 
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expression for the wrapping correction which is [50] 



e 



—2 J arcsinh 




(3) 



For the Konishi J = 2, so this gives a factor of g^. Together with the so- 
called string frame phase factors P3] which promote the string length (light 
cone size of the worldsheet cylinder) to the 'spin chain length' these give 
another factor of g'^. Together the two factors give indeed g^, making this 
contribution to appear at 4- loop. We expect therefore that the correction 
coming from two virtual particles should appear at least at order g^"^ (or if 
the string frame phase factors will also appear systematically this might be 
even as late as g^^). Therefore it is expected that the 5-loop wrapping part of 
the Konishi anomalous dimension will also appear from the leading Liischer 
correction. What makes this computation interesting, and what is the main 
motivation for our study, is the fact that two new features which were absent 
in the 4-loop case make their appearance here. 

Firstly, the dressing factor of the S-matrix between the mirror particle 
and the physical particles being the constituents of the Konishi state behaves 
like explig"^ phase), where the phase involves contribution from an infinite set 
of BES coefficients. This is in stark contrast to the behaviour of the dressing 
phase between physical particles where it behaves like exp{ig^ phase'), and 
at higher orders the higher BES coefficients enter only one by one. We will 
explicitly compute the dressing phase between mirror and physical particles 
in section m 

Secondly, the virtual particle will also modify ABA quantization condi- 
tion. The reason that this effect appears only at 5-loop order is that the 
momenta of the constituent particles get shifted from the ABA value Paba 
by a term of order g^: 



and hence their contribution to the energy coming from the dispersion rela- 
tion would appear only starting from order g^^ 



E{p) = -i / 1 + le^f^ sin^ - = . / 1 + 165(2 sin^ — — + 4 smpABA^wPABAQ^^ + ■■■ 



The appearance of these two new effects is the main motivation for our 
calculation. 



P = Paba + g ^wPaba 



(4) 




(5) 
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Of course the utility of this computation as a testing ground of the above 
two phenomena depends on the abihty of having a cross-check. For the 
Konishi operator, a direct perturbative 5-loop computation seems to be be- 
yond reach, however there are some very stringent internal crosschecks of the 
eventual final formula. 

Firstly, from the structure of perturbative gauge theory integrals we ex- 
pect the final answer to have a rather simple transcendentality structure - 
a linear combination of zeta functions (and possibly their products). Yet 
generically the subexpressions which appear when performing the computa- 
tion from the string theory side are much more complicated like polygamma 
functions evaluated at irrational complex arguments. A nontrivial cross check 
will be the cancellation of these terms between the various parts of the com- 
putation, in particular between the contribution of the dressing phase, the 
modification of Asymptotic Bethe Ansatz and the direct F-term integral. 

Secondly, purely from the integrable quantum field theory point of view 
we do not expect to have a contribution of so-called /i-terms due to the 
fact that at weak coupling bound-states are much heavier than fundamen- 
tal particles. Typically /x-terms arise from certain 'dynamical' poles of the 
F-term integrand. A consequence of the vanishing of //-term contribution 
would be a cancellation of the residues of those terms when summed over 
all bound-states. This again necessitates a subtle cancellation between the 
residues coming from the dressing phase, ABA modification and the F-term 
integrand. Therefore we may use these two consistency checks as a test both 
of our formalism and of a huge part of the BES dressing phase. Another 
motivation is the ongoing search and proposals for the nonlinear integral 
equations/functional equations which would exactly describe the spectrum 
for any size of the cylinder. The current computation could then be used as 
a test of these proposals especially as the modification of ABA quantization 
appears in a nontrivial way and the dressing phase expression in this regime 
is quite complicated. 

We will also consider a single impurity state with momentum p = it which 
can be considered as an analytical continuation of the twist-two operators 
considered at 4- loop in [32] • It should also coincide with a physical state in 
the /3-deformed theory at = 1/2. Such states have been considered pertur- 
batively in [51] and it may be possible to have a direct 5-loop perturbative 
computation in this case. Here the modification of ABA quantization is ab- 
sent, however there is a contribution from the infinite set of coefficients of 
the BES dressing phase. Therefore a direct perturbative computation would 
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be very interesting even for the single impurity case. 

3 Multiparticle Liischer formulas and ABA 
modification 

In this section we will review the formalism of multiparticle Liischer correc- 
tions introduced in [39]. In order to avoid ambiguities associated with the 
fact that the S-matrix for the AdS^ x superstring is a nontrivial function 
of both momenta and does not depend just on the difference of rapidities as is 
the case for relativistic integrable quantum field theories, we will consider the 
construction of the Thermodynamic Bethe Ansatz for a theory with diagonal 
scattering but with an S-matrix without the difference property. Then we 
will consider the construction of excited state TBA by analytical continua- 
tion along the lines of ^0]. We will then extract the appropriate formulas for 
multiparticle Liischer corrections by a large volume expansion recovering the 
expressions proposed in [39]. However we hope that the present derivation 
will make their origin clearer. We would like also to keep throughout the 
computation the conventions for the S-matrix predominantly used for the 
AdS^ X case which are different from the ones used usually for relativistic 
integrable QFT's. 

3.1 TBA for diagonal scattering 

The starting point for the derivation of the Thermodynamic Bethe Ansatz 
is the consideration of the mirror theory with space and time interchanged. 
We denote the momenta of the mirror theory by p to make a clear distinction 
compared to the momenta of the original sigma model what we denote by p. 
In comparing with the results of [30j we note that p was denoted there by 
Ptba- We consider a theory with one species of particles scattering with the 
S-matrix S{pi,p2). 

In this theory the ABA takes the form|§ 

e^P.R= J] S{p„p,) (6) 

As is standard in TBA we are interested in the thermodynamic limit R ^ oo 
of the free energy, where the physical size L is identified with the inverse 

^We use conventions used in AdS/CFT literature. 
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temperature. We will parametrize the momenta by z (which may be identified 
with the coordinate on the rapidity torus). In the thermodynamic limit, the 
roots become very dense and can be substituted by continuous densities of 
roots p(z) and densities of holes ph{z^ (unoccupied roots cf. [9]). 

Taking the logarithm of ([6]), and the derivative w.r.t. z we obtain a 
relation between p{z) and Ph{z): 

2n{p{z) + Phiz)) = Rp\z) - -d, logSiPiz), ■) * p = Rp\z) - * p (7) 

t 

In order to get a second equation necessary for solving for both densities, we 
have to minimize the free energy which is given by 

- RLf = -LH + S = -L J E{z)p{z)dz + S[p, ph] (8) 

where E is the mirror energy, while the entropy is 

S[p,Ph]=j dz{{p + ph)log{p + Ph) - plogp- phlogph} (9) 

It is convenient to introduce the pseudoenergy e{z) which is related to the 
densities of roots through 

(10) 



P + Ph I + 6-" 

and exremizing the free energy after taking into account the relation 

Sp + Sph = -^(p* Sp (11) 
/vr 

coming from ABA, we obtain the final form of the TBA equation 

e{z) = LE{z) <P{w, z) log (l + e'^^'")) (12) 

We can now evaluate the thermal free energy of the mirror theory which gives 
the physical ground state energy of the original theory at zero temperature 
but at finite size. 

E = Lf = - l^p'iz) log (1 + e-(^)) (13) 
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3.2 TBA for excited states by analytical continuation 

Initially, the Thermodynamic Bethe Ansatz could give information only on 
the exact ground state energy, and it was not known how to extract similar 
information for excited states in finite volume. In [10], a version of the Ther- 
modynamic Bethe Ansatz for excited states was constructed by making an 
analytical continuation of the ground state equations. The idea was that in 
deforming contours of integration in the TBA equations, zeroes of 1 + e"^*-^-* 
would generate poles in the integrand which could then be rewritten as ad- 
ditional source terms of the equations. Here we will perform an analogous 
heuristic procedure, fixing the sign of the pole contribution (which depends 
on the orientation of the contour w.r.t. the singularities) in order to re- 
produce ABA results at large volume. Once this is fixed, we will be able 
to unambiguously define leading corrections. Note that we assume here that 
there are no /x-terms so the calculation is analogous to the Sinh-Gordon model 
in the relativistic case and so we will assume that each physical particle is 
represented just by a single pole. 

We will need the relation between the mirror energies and momenta {E, p) 
and the physical ones {E,p). These are defined by 

E = ip p = iE (14) 

which is fixed by taking the mirror integral to be for z = x — 002/2 with x 
real, and consequently identifying z^ / z'^ with e~^. An opposite choice is 
also possible. 

Let us first consider the energy formula ([13!) integrated by parts: 

E = I ^p{z)dJog{l + e^^^^^) (15) 

Suppose that 1 + e"^*^^!'^) = g. Then the orientation has to be such that the 
dz log (1 + e"^(^)) contributes —1 by residues. We thus obtain 

E = E{z,) + E{z,) - j |^p'(;,)log(l + e-(^)) (16) 

By the same mechanism an analogous contribution will arise from the TBA 
equation (fT2|) . with the sign already fixed by the above considerations. We 

get 

^(9^1og5(w;,z))log(l + e-^("')) 

(17) 
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The multiparticle Liischer formulas will follow, for this theory, by performing 
large volume (large L) expansion. 



3.3 Multiparticle Liischer formulas 

Firstly in order to compute the energy keeping the first exponential correc- 
tions, we may neglect the integral term of (fT71) when inserting e{z) into (fT6|) . 
We thus get 

E = E{z,) + E{z2)- I ^p'e-^^ ^ 



27r S{zi,z)S{z2,z) 
= E{z,) + E{z,) - 1 ^ e-^^5(^, z,)S{z, z,) (18) 

We recognize at once the F-term integral (with q = p). 

In order to complete the formula we need to self consistently fix the 
positions of the poles Zi and Z2. We will do it in two steps. First we neglect 
the integral term of f|T7|) and impose for the rest 6{zi) = in + (27rn)i. We get 

in = e{zi) = iLpi + in + log5'(2;2, Zi) (19) 

where we supposed that S{z,z) = —1. This gives at once ABA equation 
with our conventions 

e^^P^ = 5(pi,p2) (20) 

However it turns out that we have to be more precise in the determination 
of the position of the roots and we have to include the integral term. In order 
to define the quantization conditions we thus have to use 

e{z) = iLp{z) + log S{zi,z) + log S{z2,z) + 

(21) 

J 2m S[w,z) 
The quantization conditions s{zi) = in takes the form 

= log{e''^P^S{z2,z,)}+ [ p^{d^Siw,z^))S{w,Z2)e-'^^^-^ (22) 
^ ' J 2m 

BYi ' V ' 

$1 

= log{e''^P^S{z,,Z2)}+ [ p:S{w,z,){d^S{w,Z2))e''^^'~'"^ (23) 

^ ' J 2m 

BY2 ' V ' 

$2 
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Since the integrals are exponentially small we may solve 
terms of corrections to ABA giving 



these 



equations in 



dBYi 

dpi 
dBY2 
dpi 



Spi + 



Spi + 



dBYi 

dp2 
dBY2 

dp2 



Sp2 + $2 = 



Sp2 + $ 







(25) 



(24) 



The final formula for the energy thus takes the form 



E ^ E(pi)+E(p2)+E'{pi)5pi+E'{p2)5p2- / e-^^5(z, ^i)5(^, ^2) (26) 



For the case of the AdS^ x string sigma model wc just have to replace the 
product of the S-matrices by an appropriate supertrace. Thus the integrand 
becomes essentially the transfer matrix. However the part coming from the 
modification of ABA quantization does not have a-priori such a direct rela- 
tion to the transfer matrix. Hence in the following we will evaluate it directly 
from its definition using the S-matrices and their derivatives. 



4 The dressing phase in the Luscher kinemat- 



As we explained in section 2 one of the main sources of corrections which 
appears at five loops is the contribution of the dressing phase. In this section 
we present two calculations for its leading order part in the Liischer kine- 
matics. By this we mean a kinematics which is relevant in calculating the 
finite size correction, that is when the first argument is in the mirror region 
|xi| < 1 while the second is in the physical one \x2\ > 1. 

The dressing phase when both arguments are in the physical region, 
(1x^2 1 > 1) can be written in the following form 



0{xi,X2) = xi^i, ^t) - Xi^i, ^2 ) - x{.Xi , ^t) + Xi^i , ^2 ) (27) 



ICS 



where 



00 



Cr,s{,g) \ 1 



1 



x(xi,x2) ^-^Y^ 



(r-l)(5-l) [x\-'xl- 




(28) 



r=2 s>r 
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The coefficients Cr^sig) have a convergent weak couphng expansion 

= 2co= {^is -r-l))(r- l)(s - 1) f ,t JMJ.^2gt) 

As was shown in [29j the dressing phase has also an alternative double integral 
representation, which is suitable for analytical continuations: 

XiXi,X2) =i(f ^ l2{Wi,X2) (30) 

Jc^ 2m wi - xi 

where 

r dw2 1 T{l + ig{wi + w];^ -W2-W2^)) 

l2{wi,X2)=f — log— — ■ 31 

Jc-i 27r« W2 - X2 r(l - tg{wi + - - )) 

and the integrations go over the unit circles. In the next subsections we 
focus on the weakly coupled regime, review the dressing phase for physical 
particles in this regime and calculate the analytical continuation for both 
representations in the Liischer kinematics. 

4.1 BES dressing phase for physical particles 

In the g ^ limit we can expand the Bessel function as 

U2gt) = ^-^{l-^ + ...) (32) 
nl n + 1 

and perform the integration in (!29l) . This provides the leading order be- 
haviour of Cr^sig)'- 

= ' r2)?^ ^"^"' + ^ - '^'^^^ + ^ - + ^^'^ 

- — (r + s)(r + s - l)!C(r + s) + . . . 1 
rs 

In calculating the Liischer correction we have to analyze the dressing phase 
(127|) for the case when the parameters Xi,X2 in the g ^ limit behave as 

q + iQ . 2g ^, , 2u±i 

2g [q - iQ) 2g 

(34) 
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Let us start with the case when both parameters are in the physical 
kinematics Xi^2 = 0'i,2g~^ + 0{g). Then the leading order behaviour comes 
from the term 02,3 ((7) and is given by 

Xix„x2) = -i2g\{3) + 20A(5))^^ + 0{g'') (35) 

The next coefficients 02,5(5') and 03,4(5) only enter at 5^° order. For physical 
particles, the dressing phase starts at order and this is what we have to 
use in the calculation of the ABA. 

Let us now describe the regime where one of the particles is in the Liischer 
('mirror') kinematic which is relevant for the five loop wrapping computation. 



4.2 BES dressing phase in the Liischer kinematics 

Suppose now that the ffist argument is in the mirror kinematics Xi = a^^g + 
0{g^) while we keep the second in the physical one X2 = a2g~^ + 0{g). Since 
Cr,s{g) scales as g^~^^~'^ while * ^ the leading order contribution 

comes from the r = 2 part of the ffist sum: 

X{x^,x2) = Y: + 0(9') (36) 

^ (S - 1) Xi X2 

Clearly this is in stark contrast to the case of the ABA as in this Liischer 
kinematics infinite number of coefficients contribute. Using the leading order 
term of the explicit weak coupling expansion of C2,s determined in (|33il we 
obtain 

Xix„x2) = ^^2cos(^(.-3))ariC(^)+0(5^) 

^ s>2 

= 9^^s,{-^a^) + S,i^a^)) + 0ig') 
a2 

where Si{n) = ^22=1 1 harmonic number which has an analytical 

continuation in terms of the digamma function i{j{x) = as 

X{xi,x2) = — (27^ + ^(1 -zai)+^(l + zai)) + 0(5^) 
02 

= ^{2^E + ^Hai)+4j{ta,)) + 0{g') (37) 

^2 
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where 7^; is the Euler constant. Thus for sufficiently small Oi the original 
BES dressing phase provided a convergent expansion for the leading order 
(7 — > behaviour in the Liischer kinematics. The resulting expression then 
can be analytically continued for any ai. Let us check that this is the right 
analytical continuation by calculating the weak coupling expansion of the 
analytically continued dressing phase from the DHM integral representation. 

4.3 BES dressing phase in the Liischer kinematics from 
the DHM integral formula 

The proper analytical continuation of the dressing phase fl30l) into the Liischer 
kinematics, compatible with crossing symmetry, was determined in 



xixi,X2) = i(f ^ hiwi, X2) - ihixi, X2) (38) 

Jc^ 2m wi - xi 

As in the previous subsection we are interested in the weak coupling ex- 
pansion of X ill the Liischer kinematics so we take Xi = a^^g + 0{g^) and 
X2 = Ci2g~^ + 0{g). Observe that since \x2\ > 1 the contour never encircles 
X2- Let us focus on the second term. We change the integration variable to 
U2 = gw2 and write 

r du2 I T{l + ig{xi + xl^)-iu2-ig'^U2^)) , . 

h{xi,X2)=f — log— — _i , . 2 -UN (39) 

2mu2 - a2 T{1 - ig{xi + x-^ ) + iU2 + ig^U2 )) 

Since the integration for U2 goes over a shrinking circle of radius g —>■ we 
pick up the contributions of the poles at M2 = only (not at 02). From the 
expansion of the F function in g^U2^ one can observe that higher order poles 
contribute at higher orders in g"^. Thus when we focus on the leading order 
behaviour we are allowed to keep the first order pole: 

J2(xi, X2) = i (^^^ [V^(l + ig{xi + x^')) + V^(l - igixi + x^'))] + 0(^^) 

(40) 

Taking into account that Xi = a-^ ^g + 0{g'^), the leading order part turns out 
to be ^ 

hixi, X2) = t mi + zai) + ^(1 - ta,)] + 0{g') (41) 
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In addition to this term we also have a term of the form 



Il{Xi,X2) = i i 



dwi 1 



l2{Wi,X2) 



(42) 



27ri Wi — Xi 



Since Xi = a^^g + 0{g^), it hes inside the integration contour. We need 
a convergent small coupling expansion. For this we change the integration 
variable to ui = wig~^ and close the contour around oo. The leading oder 
result comes only from the pole at oo and reads as: 



In order to determine the leading order behaviour of x we have to combine 
h and h- The result provides the leading order asymptotics of the dressing 
phase in the Liischer kinematics 



which is consistent with what we obtained from the BES dressing phase. 
This is the main result of this section. We will need to take into account 
this expression when we calcTilatc the five loop anomalous dimension of the 
Konishi operator in the next section. 

5 The Konishi computation 

In this section, which is the main part of the paper we describe and evaluate 
the various ingredients which together contribute to the 5-loop wrapping cor- 
rection. For completeness wc first describe the evaluation of the Asymptotic 
Bethe Ansatz (ABA) contribution to the anomalous dimension and then pro- 
ceed to compute the various parts of the wrapping correction, summarizing 
the previous 4-loop result in the process. 

5.1 Asymptotic Bethe Ansatz for the Konishi operator 

In this subsection we calculate the scaling dimension of the Konishi operator 
from the ABA. In principle the complete scaling dimension has the form 




(43) 




(44) 



A(5) = AABAig) + 



(45) 
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where the first term contains the result of the ABA, while the second the 
wrapping contributions. In doing the calculation we choose the representative 
for the Konishi operator in the s/2 sector. ABA amounts to calculate in this 
case the energy of the two particle state in volume L = 2 with momentum p 
and —p which satisfy the equation 



e 



ipL 



S{p, -p) = _pyie(xip)M-P)) (46) 

where S^''^'^\p, —p) is the scattering matrix in the s/2 sector which reads as 

and x^{p) = ^(cot(|) ± + E{p)). We can solve this equation perturba- 

tively. We take the Ansatz p = '^^p^^^'g'^^ and systematically expand both 
the scalar and the dressing part in increasing orders of g. The first two orders 
can be determined without the dressing phase. These solutions then can be 
put into formula ( |35ll to determine the leading contributions of the dressing 
phase. Finally we find the solution for p at the required order 

P = f-V3g'+^-^g'-24V3{l + Cm9' (48) 

+ ^(671 + 960(C(3) + C(5)))/ + Oig'') 
This momentum then determines the ABA part of the scaling dimension as 

AABAig) = 2E(p) =4+12/-48/ + 336/-(2820 + 288C(3))/ 

+ (26508 + 4320C(3) + 2880C(5))^^° + 0{g^^) 

The wrapping part of the anomalous dimension starts as 

A.(^) = ALV + ArV° + 0(^^^) (49) 

In the 4-loop part has been evaluated to give 

A^^'^ = 324 + 864C(3) - 1440C(5) (50) 

and the aim of the rest of this paper is to calculate the wrapping part at 
5- loop order Aw^\ Our final result is 

A(J°) = -11340 + 2592C(3) - 5184C(3)^ - 11520C(5) + 30240C(7) (51) 
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5.2 The general structure of the wrapping correction 



There are two sources of the Liischer correction for the Konishi operator 
which have to be taken into account when we calculate the 5-loop anomalous 
dimension (recall formula (l26l) ). 



Modification of Asymptotic Bethe Ansatz quantization 

First of all the leading contribution which comes from the modification of 
ABA is of order g^^. This has the form 

At'''' = E'ip^)6pi + E'ip2)Sp2 (52) 

where pi = —p2 = p is the solution ( HHl) of the ABA, while 6pi are the shifts 
due to the virtual corrections. For the 5-loop result it will be enough to take 
p = 27r/3 in this formula. 

The shifts 6pi for the Konishi operator can be found from relations (12^ 
and fl25p which in our case take the form 

- yp2 + $1 = 
-^Spi + —5p2 + ^2 = 
where $ = $i = —$2 and $ for the Konishi operator is as follows 

Solving the equations for 6pi and 6p2 we found that 

Spi = -Sp2 = 

This means that the momenta of the Konishi constituents are shifted in 
opposite directions by the same factor. It leads to the vanishing of the total 
momentum after ABA modification as it should be. 
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When plugged into the dispersion relation it gives us the energy shift (1521) 
due to the Asymptotic Bethe Ansatz modification as 



At^^ = -1$ (54) 
We will evaluate $ explicitly later in the section. 
The F-term integral 

The second contribution to the 5-loop result comes from the expansion of 
formula [39] 

up to subleading terms. We parametrize the virtual particle using the mirror 
momentum q as in [39] (see also fl68l) ). while for the physical particles we 
use rapidities as in [3.2j. Hence Ui = —uu = u are the rapidities of the two 
particles forming the Konishi state. They are related to the momentum p as 



The rapidity variable for the Konishi state is therefore given as 

M = Mo + U2g'^ + ... = ^ + -^^^ + . . . (57) 



To simplify our further considerations we rewrite ( !55l) in the form 

OO POO 7 - 

= -E/ ^{^) Soiq,Q,u)S^iq,Q,u)Sdq,Q,uy (58) 

where we split the contribution coming from the S-matrix into three parts: 
the scalar, the dressing and the matrix part. Each of those parts involves a 
product over contributions of the two particles forming the Konishi state. 
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It is convenient to perform the weak coupling expansion of this formula 
in two steps. First we keep variable of order 1 and expand all other 

sources of g dependence. Then we take into account the fact that u itself is 
g dependent and compute the g dependent terms due to fl57j) . 

All ingredients in formula (l58l) can be expanded into a series in g"^, keeping 
u fixed, as 



z 



So{q,Q,u) = Sy{q,Q,u)+g Sy{q,Q,u) + ... 
S,{q,Q,u) = l+g^S!^\q,Q,u) + ... 
Sm{q,Q,u) = g^S^\q,Q,uf + g^S^\q,Q,u) + . . . 

Now we will factor out the leading piece 

Y^'\q,u)g' = s!,'\q,Q,u)Sg\q,Q,ufT('\q,Q,u)g' (59) 

and rewrite the first subleading term as a sum of contributions coming from 
the matrix part, scalar part, the exponential term and the dressing factor: 



2 m \ -J i -V ; - / _|_ 



In the following we will calculate all these contributions one by one. 
Since YQ^\q, u) is already at order g^^, we may safely set here u = 

Y^''\q^u)g'^ = Yi^'^'\q)g^^ + 0{g'^) (60) 

where y^^°-°)(g) = y^^°)(g,«o). 

As mentioned earlier we also have to take into account that the rapidity 
u is g dependent by itself. Clearly the only part that will contribute at order 
g^^ is the expansion of fl59|) 



Y^'\q,u)g^ = yi^'°)(g)/ + Yl,'''\q)g'' + . . . (61) 
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with Yq'^^ (q) = Yq^ (g, Uq) ■ To summarize, the wrapping correction at 4- 
and 5-loop respectively take the form 



=1 

oo 



= - E £ S + Yr\^) + yt\^)) (63) 

The leading contribution and 4-loop wrapping 

The leading expansion of integrand YQ{q,u) is known to be proportional to 
and was found in [39j as 



Y^'\q,u)g' = sS'\q,Q,u)S^^\q,Q,uyT^'\q,Q,u)g' 
lQ384g^Q^{-l + q^ + - Au^ 

1 



((g - 2(g - 1))2 - Au^){{q - i{Q + 1))2 - Au^) 



When expanded further using the fact that the rapidity variable u is 
(yf^-dependent ( 157|1 it can be rewritten as 

Y^'\q,u) = Y^'''\q)+g'Y^'''\q) (64) 

where Yg^^'^^^q) is the only contribution which is relevant for the 4-loop 
calculations. 

The integrand is a rational function and hence the integral can be carried 
out by residues. The poles are of two kinds. Firstly there is a fourth order 
pole at q = iQ which we call 'kinematical' pole since it comes just from the 
exponential terms. The remaining poles come from the S-matrix and are 
'dynamical' i.e. s and t channel poles. They would be associated to possible 
/i-terms. Since, as argued in [39], /z-terms should not be present at weak 
coupling, we expect that the residues of the dynamical poles should sum up 
to zero after summation over all Q. We can indeed verify that this is the case. 
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Hence the 4-loop wrapping correction to the Konishi anomalous dimensions 
can be obtained just by summing the residue of the kinematical pole only: 

oo 

A® = -27ri V res Y^^'°\q) = 324 + 864C(3) - 1440C(5) (66) 

^ — ' q=tQ ^ 

Q=l 

We expect that the same cancellation of 'dynamical' poles will also hold at 
5-loop. Indeed this will be a nontrivial internal consistency check of our 
calculation. 

5.3 The 5-loop integrand 

The full integrand which will be relevant for 5-loop calculations consists of the 
subleading contribution of YQ\q,u) (namely YQ''^\q)), the leading contri- 
bution of YQ^°^(g, m) (namely YQ^''^\q) which comes from the matrix, scalar, 
exponential part and from the dressing phase) and the leading contribution 
^giq) coming from the ABA modification. We will now evaluate one by one 
the individual contributions. 

Matrix part 

The subleading matrix part of integrand S^\q, Q, u) for the Konishi operator 
can be evaluated using the formulas (78)-(82) from [39j. When normalized 
with respect to leading term {q,Q,u) it stands as 

S^i\q,Q,u) 16g(l -iq-Q- Au'') 16 

S^i\q,Q,u) ~ {q + ^Qm-^{Q-l)f-Au^){l + Au^) (l + 4nT 
16(g^ + g^ - 1) 4(-5 + + 4Q + 5Q^ - 4^^) 

Scalar part 

The scalar part of the integrand can be evaluated using the S-matrices of the 
scattering of the si (2) bound-state constituents (^rf , z^), . . . , (zg, Zq) with 
the fundamental magnon (a;~, x^): 

Q 

Soiz^,x^) = llS^^^'\zt,x^) (67) 

i=l 
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with 



Z+ — X 1 



In order to calculate (1671) we use for z-^ the value z and for Zq the value 
z~^ where 

(68) 




whereas z^ we take from 




The intermediate are determined from the pole condition 2;^ = z^_-^^. In 
principle we could choose different signs in front of the square roots. This 
amounts to choosing a different representative for the constituents of the 
Q bound-state. It was shown however in [12] that all of the different con- 
stituents lead to the same S-matrix after analytical continuation. Our choice 
is technically the simplest one and leads to the expansion of the bound-state 
parameters up to the second order as 

^+ ^ 2tk + q-tQ _2g 2g 2g 

^ 2g q-iQ 2ik + q-iQ q + iQ ^ ^ 

The result for the leading order scalar part of the integrand is then 



{{q - t{Q + 1))2 - 4u^){{q + t{Q + 1))2 - Au^) 

1 



((g + z(g-l))2-4n2) 
while the subleading one can be split into two parts: a rational part 

Slil{q,Q,u) ^ 32 32 

Sjl'\q,Q,u) {^q + Q){l + iu') ^ (1 + 4«T 

8 / 2g(g-z(g-l)) _ 2g(g + z(g - 1)) 
g2 ^Q2{^q_ _ i))2 _ 4^2 + _ i))2 _ 4^2 + 

2q{q-t{Q + l)) _ 2q{q + t{Q + 1)) _ 2{q^ + Q"" + 2Q) \ 
(g-i(g + 1))2 -4m2 (g + ^(g + i))2 _4y2 1 + 4m2 ) 
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and one which contains polygamma functions 
Exponential part 

The leading and subleading term of the exponent which appears in Yq is 
found to be 



where and z~ is taken from 



1-9 



2 16 



(g2 + g2 



Dressing part 

The dressing part can be found from formula fl37|) as 



S!^\q, Q, u) = -Y^(7i. + ^^(^(-^? - Q)) + l^ili^Q + Q))) (71) 
ABA modification 

To make our notation more compact we rewrite the leading order Asymptotic 
Bethe Ansatz modification formula into the form 

^(8) = r ^ u) = Y,r ^Y^'\q, n)ABA(g, Q, u) (72) 

In order to find ABA{q,Q,u) it is enough to use formulas (78)-(82) from 
[39j . Following fl53l) we have to take the derivative of S-matrix elements with 
respect to q and then calculate the supertrace. The final result we obtained 
is 

ABA{q,Q,u) = — ^^^ + — ^-7m^+ ^ 



+ 



q^ + Q"^ —i — q — iQ + 2u i — q — iQ + 2u 
1 1 2(g + 2u) 



-i-q + iQ + 2u i-q + iQ + 2u I - q^ - Q'^ + Au^ 
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The contribution YQ''^\q) 

As mentioned in tlie previous section Yq ' (g) comes from tlie 4- loop inte- 
grand expressed in terms of u, when we take into account the shifts of the 
rapidities of the constituents of the Konishi due to ABA. In order to find it 
we have to plug u = + ■^g'^ into YQ'\q) and expand it to the second 
order in g"^. The result is 

(8,2). ^ ^ 4718592g^(16 + V - 12Q^ + 9Q^ + 6g^(-10 + 3Q^)) 
^3 (9g4 + g^2(2_6Q + 3Q2) + (4-6Q + 3Q2)2)2 ^ 

(-4 + 3g2 + 3g2)(_i6 _ 9g4 _ i2Q^ + 27Q^ + Qq\-2 + ?,Q^)) 



(9g4 + 6g2(2 + + 3g2) + (4 + eg + 3g2)2)2(g2 + g2)4 
5.4 Integration 

Before we proceed, it is fruitful to observe that we can symmetrize the inte- 
grand with respect to q without changing the result of integration 

When we calculated the integrals we symmetrized some parts of the integrand 
leaving the rest not symmetrized depending on which form is easier to handle 
and gives simpler result. 

Apart from the polygamma functions which appear in the dressing and 
scalar part, the remaining part of the integrand is a rational function. It 
can be then integrated over the real line by taking residues at the position 
of poles lying above the real line. All such poles can be classified into two 
groups: poles coming from the S-matrix parts of the integrand ('dynamical' 
poles), four of which lie above the real line: 

g = ^(g±i)±i=, (74) 

and a pole coming from the exponential part which is exactly q = iQ. 

It turns out that when we symmetrize the whole integrand the contri- 
bution coming from dynamical poles vanishes when summed over Q as in 
the 4-loop cas It can be explained by the fact that at weak coupling we 



-"^"It can be checked both numerically and completely algebraically using the methods of 
the subsequent section. 
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expect the //-term contributions to be absent. We want to stress that this is 
a far from trivial consistency check of our formulas. 

Beside the residues of the rational function poles we have to handle the 
additional poles coming from the polygamma functions which appear in the 
dressing and scalar part. In order to simplify our considerations, it is conve- 
nient to symmetrize the S^(q, Q, u) part. It turns out that in that case the 
polygamma functions will vanish leaving us only with rational function 

(^(1 -«)- >*(iH, + «)))'"■" 

So the only remaining polygamma functions come from the dressing part. 
The positions of polygamma function poles are known to be at negative 
integers. In our case it means that j{—iq — Q) = —n or ^{iq + Q) = —n 
where n > 0. Solving with respect to q we obtain the positions of integrand 
poles coming from polygamma function as 

q^i{Q-2n) n > or q^ i{Q + 2n) n>0 (75) 

We want to close the contour of integration over the real line meaning we have 
to take into account only those residues which have nonnegative imaginary 
part. 

In the proceeding sections we will calculate all the residues. We have 
chosen to symmetrize only two parts of integrand: the scalar integrand part 
which contains polygamma functions and the part coming from the ABA 
modification. The rest is left in a nonsymmetrized form. 



Residues of rational functions 

For the rational part of the integrand the result after taking the residues 
(regardless if we take q = iQ or 'dynamical' poles) can be rewritten as a sum 
of terms with minimal denominators which are of the form 

which will give zeta functions when summed up or 

aQ + b 



1 ± 3Q + 3g2 
which will produce polygamma functions. 



(77) 
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For the dressing part the residue at q = iQ gives exactly 



3456(-2 + 12Q^ - 45Q^ + 27g'^)C(3) 

Q3(l_3g2 + 9Q4)2 ^ ) 

while the residues at the dynamical poles are some rational functions of Q 
with complicated coefficients containing '?/'(^(3 + i^/S)) and ^/'(|(3 — 



Residues of polygamma functions 

To calculate the contribution coming from the polygamma functions appear- 
ing in the dressing part of the integrand we want to find all the poles of the 
form ( l75l) which have nonnegative imaginary part. It is easy to notice that 
such poles are of the form gQ^„ = i{Q + 2n) where n is (possibly negative) 
integer obeying Q + 2n > 0. The residues at the positions qq^n is found to 
be 

-864g2(i + 3r^2 + 3nQ)hign{n) 
n^{l - 3n2 + 9n^){n + Q)^{{1 + 3n^ + 6nQ + 3Q^y - {3n + 3Qy) 

The remaining task is to sum the above formula over n and Q. We have to be 
careful during the summation process because there exist one pole which lies 
on the real line for every even Q. In that case we have to take only half of the 
residue with the minus sign due to how we have chosen the orientation of our 
integration contour. Additionally, we have to remember that the residues at 
q = iQ were taken into account before. Keeping it in mind we calculated the 
sum over Q and obtained 

1728(3^2 - 1)(-1 + 3^2 + 3n^ - + 9n^ + {n^ - 3n^ + 9n'^)^^^\l + n)) 

^6(1 - 3^2 + 9n4)2 

which have to be summed over n from 1 to oo. It can be done using methods 
from the next section. 



6 Summation over bound-states 

There are two types of sums over bound-states appearing in our calculations 
after the residues are taken. Firstly, there are sums of the form fl76l) or flTTl) 
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which can be easily summed using Mathematica. The first ones give us (- 
functions while the second ones contain polygamma functions of the form 
7/^(1(3 + iVS)) and V^(|(3 - ^^3)). 

On the other hand during the 5-loop Konishi calculation more difficult 
sums emerge 

oo 

E(™) = J2 R{Q) ^^""\Q) m>0, (79) 
Q=l 

where R{x) is a rational function of x and ^^"^\x) is the mth polygamma 
function given by the definition 

^M(a;) = ^!l^, m>l, (80) 

with ^{j(x) = ip^'^\x) being the digamma function i{j{x) = ^ J^^^ ■ The 
evaluation of the sum fl79l) goes as follows: R{Q) is decomposed as a sum of 
two terms R{Q) = Ro{Q) + Ri{Q), where Ro{Q) contains the sum of pure 
power terms of the partial fraction decomposition of R{Q) (i.e. Ro{Q) = 
+ ^2 + •■O' while Ri{Q) contains the rest. In this case the sum fl7^ is 
decomposed into two parts as well: 

^M^^^M^^M^ (81) 

where 

oo 

Si-) = ^i?,(Q)^('")(g) m>0, a = 0,1. (82) 

Q=l 

Using a series representation for the polygamma functions, Sq"*^ can be ex- 
pressed in terms of the values at infinity of nested harmonic sums, which 
can be expressed in terms of multivariate zeta functions (Zagier-Euler sums) 
|43j . These sums can be reexpressed in terms of ordinary Euler sums. The 
relations between the former and the latter can be found using EZ-Face - 
an online calculator for Euler sums 01]. At the end of the process Sq™-* is 



^0 

expressed in terms of zeta functions taken at integer values. 

The calculation of H^i^^ requires a different method: representing Ri{x) as 
the Laplace transform of its inverse Laplace transform, and using appropriate 
integral representation for the polygamma functions, S^™'' can be transformed 
into a double integral form, in which the summation can be easily performed 
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and the remaining double integral expression can be calculated exactly with 
the help of Mathematica. For details see Appendix B. 

Using the methods presented above all the sums can be evaluated. The 
striking observation is that all 'nasty' polygamma functions disappear leaving 
us only with the (^-functions. The final result is 

A(]°) = -11340 + 2592C(3) - 5184C(3)2 - 11520C(5) + 30240C(7) (83) 

As for the 4-loop case the result is a sum of (^-functions of odd degrees with 
integer coefficients and correct transcendentality degree. The new feature of 
the result is that products of ^-functions start to appear which have been 
absent in the 4-loop case. 

To summarize the anomalous dimension of the Konishi operator up to 
5-loops is 

A = 4+12/-48^^ + 336/ + 96(-26 + 6C(3)-15C(5))/ (84) 
-96(-158 - 72C(3) + 54C(3)2 + 90C(5) - 315C(7))^^° 

7 Conclusions 

Anomalous dimensions of operators in A/" = 4 SYM correspond to energies of 
string states in AdS^ x 5*^. The leading perturbative orders are given by the 
Asymptotic Bethe Ansatz, while at a certain loop order, there appear new 
contributions coming from a topologically distinct class of Feynman diagrams 
- so-called 'wrapping interactions'. On the string theory side these corre- 
spond to virtual particles propagating around the string worldsheet cylinder. 
The leading corrections arise from multiparticle generalizations of the classi- 
cal Liischer terms. Further corrections are due to many virtual particles and 
all these should in principle be resummed by a TBA system. 

The most convenient testing ground for these issues is the shortest non- 
protected operator in A/" = 4 SYM - the Konishi operator, as well as single 
impurity operators in the /9-deformed theory. 

In a previous paper |39], the leading four loop part of the wrapping cor- 
rection has been found. It came from a two-particle Liischer term and gave 

^^KomsM = 324 + 864 C(3) - 1440 C(5) (85) 
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The aim of the present paper was to compute the five loop wrapping con- 
tribution. Although it still arises just from the two-particle Liischer term, 
it is very interesting as it involves two new ingredients. Firstly, there is a 
nontrivial modification of the Asymptotic Bethe Ansatz quantization due to 
the sea of virtual particles. This term is quite interesting as it does not seem 
to be simply related to the transfer matrix appearing in a direct expansion of 
the proposed TBA systems. Secondly, the BES/BHL dressing factor cannot 
be neglected any more, and moreover, due to the specific kinematics of the 
Liischer term, an infinite set of coefficients of the dressing phase starts to 
contribute at once. The result obtained in the present paper is 

<ln../.. = -11340 + 2592 C(3) - 5184 C(3)' - 11520 C(5) + 30240 C(7) (86) 

This computation is subject to two nontrivial cross-checks. Firstly, the par- 
tial results (coming from the dressing phase, modification of the ABA and 
the remaining part of the S-matrix) have a very complicated transcendental- 
ity structure involving polygamma functions. All these cancel out leaving the 
final result as a simple combination of odd ( functions. Secondly, the residues 
of the 'dynamical' poles (associated with /i-terms) should cancel out between 
the various terms when summed over the types of bound-states. Both of these 
cross-checks are satisfied in our case and involve a rather intricate conspiracy 
between the various terms. 

In Appendix A, we have computed the five loop wrapping correction to a 
single impurity in the /3-deformed theorj0. In this case the four loop result 
was (M = 1 in [32], see also [I5lli6] ) 

ASU;e = 496C(3)- 640C(5) (87) 

while the five loop wrapping contribution computed in Appendix A is 

<sUe = -1536 C(3)^ - 4096 C(3) - 5120 C(5) + 13440 C(7) (88) 

The ABA modification does not appear here, however the infinite set of BES 
coefficients contribute just as in the case of the Konishi operator. 

It would be very interesting to verify these results perturbatively. Espe- 
cially the single impurity operator may be within reach. This might help 

^^This computation relies on an additional assumption on performing analytical contin- 
uation from even to odd spins in [32j . The four loop result has been verified by a direct 
perturbative computation in |41j . 
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in understanding the precise relation between the structure of direct per- 
turbative computations and the string theory computations based on the 
S-matrices and Liischer formulas. 

Another application of these results would be to test the excited state 
TBA systems recently proposed for AdS/CFT. It would be very interesting 
to see the appearance of the rather intricate modification of the Asymp- 
totic Bethe Ansatz and the specific analytical continuation of the dressing 
phase from these formulations, especially as both of these ingredients are very 
strongly constrained by the transcendentality structure and the cancellation 
of dynamical poles between themselves and the rest of the integrand. The 
four loop result is sensitive to the TBA source terms and has been rederived 
in [55]. The five loop result depends on the structure of the convolution 
terms and thus is a more sensitive test of excited state TBA systems. In 
particular the procedure of section 3 to obtain the Liischer corrections could 
be applied e.g. to the TBA system in the second paper of [35]. However 
this is technically quite involved. The five loop wrapping correction for the 
Konishi operator seems therefore to be an interesting and robust test for the 
excited TBA systems. 
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A Liischer correction for a one particle state 



In this Appendix we calculate the subleading wrapping correction to the en- 
ergy of a one particle state. Such a state must have vanishing momentum 
in the Af = 4 super- Yang-Mills theory thus its energy is protected. In the 
P deformed theory, however, at /5 = | a one particle state with vanishing 
rapidity u = {p = tt) is allowed and acquires nontrivial finite size correc- 
tions. The condition for its vanishing rapidity is protected at all orders both 
in the ABA and in the wrapping corrections. So in this case we will not 
have a contribution coming from the modification of the ABA quantization. 
However there will be a nontrivial direct wrapping correction to the energy, 
which will include, starting from five loops, a contribution from the BES 
dressing phase. Thus the single impurity operator provides a testing ground 
for that part of the Liischer correction. 

The leading order wrapping correction can be written as 



oo „ , oo „ , 

Q=l-^ 0=1 



(89) 



Here SqI represents the scattering matrix of the mirror Q particle with the 
fundamental u = physical particle. We can further decompose the scatter- 
ing part as the scalar part (5*0), the dressing part (S'o-) and the matrix part 

S2l{q, 0) = So{q, Q)S^{q, Q)Sm{q, Qf (90) 
We expand each quantity to subleading order in g"^ as 

So{q.Q) = Sl,'\q,Q)+g'S^'\q,Q) + ... 

S.{q,Q) = l + g'Si'\q,Q) + ... 

Sm{q,Q) = g'S^i\q,Q)+g''S^^\q,Q) + ... 

When we calculated the leading wrapping correction to twist two operators 
for odd particle number [32], we observed that in order to be compatible with 
gauge theory calculations and to provide the proper analytical continuation 
from even cases to odd ones we had to omit the contribution of the fermions. 
We accept this convention now too, but call the attention for the need of a 
derivation of this proposal based on first principles. Under this assumption 
the leading order correction of the matrix part is 

{q' + Q'){Q + ^q-l) ^ ^ 
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while the subleading one is 

S^^\<1,Q) _ i{q' - ^q^Q - 3) - ziQ - If + q{Q' + 1)) 



Si\q,Q) {q' + Q'){q-^{Q-l)) 

The leading and subleading correction of the exponential part is given by 



(92) 



e-'^'^^^g'T^'\q,Q)+g'T^'\q,Q) 



16g' 



(g2 + g2)2 



9 



16 



(g2 + Q2) 



(93) 



In calculating the scalar part we have to use the parameterization of the 
bound-state as we did for the Konishi operator. The result for the leading 
order scalar part reads as 



SS'\q,Q) 



q-t{Q-l] 



{q + t{Q - l)){q - t{Q + l)){q + t{Q + 1)) 
while the subleading contains a rational part 



(94) 



sS'\q,Q) 



= 8 + 



16(-1 + 2iq)q'^ - 16(2 + q{i + 6q))Q 

(1+V)(^2 + Q2) 

16qii-l + 2(f + Q) 2q(3 + 2Q) 



(1 + V)(^2 + (_1 + g)2) (1 + 4^2)(^2 + (1 + g)2)) 

and a polygamma part: 



5 fo O) 1 1 

Jy)^ = 8V'(--(ig + Q)) - Si^i-i-iq + Q)) 



The dressing part reads as 



27i? + ^(-^(^5 + Q)) + V^(^(^? + Q)) 



The full integrand can be written as 

YQ{g) = gX\<l) + 9'X'\q) + . . . 
The leading order part is given by 

y^'\q) = s'^o\Q.Q)s'i\q,Qf^^'\q,Q) 

4096g2(-i + ?2^g2)2 



(95) 



(96) 



(97) 



{q^ + g2)^(?^ + (-1 + g2)2 + 2q\i + g2)) 
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It has a kinematical pole at iQ and two dynamical poles at i{Q ± 1) on 
the upper half plane. If we take the residue of the integrand only at the 
kinematical pole and sum over Q we obtain the leading wrapping correction 

A(„«) = 128(4C(3) - 5C(5)) (98) 

We note that the contributions of the dynamical poles summed over Q cancel 
out, so this is indeed the full leading order result, which has been verified by 
a direct perturbative calculation on the gauge theory side in [IT] . 
In calculating the 5-loop subleading correction we can write 



We analyze separately the contributions of the rational part and the polygamma 
part. We can observe that the rational part of the scalar contribution to- 
gether with the matrix contribution and the exponential part gives a sym- 
metric function in q: 

(8(g^(7 + 2Q) + (3 + 2g)(-l + Q^)^ + 2q\h + 2Q){1 + Q^)) 
(g2 + (-l + g)2)(g2 + g2)(g2 + (i + g)2)) 

Taking the residue at the kinematical pole and summing over Q gives 

— (727r2 - 307r^ + 27r^ - 1296C(3) - 1080C(5) -|- 2835C(7)) (99) 

The contributions of the dynamical poles do not cancel when summed over 
Q (as we did not analyze yet the full expression) and give 

_^(_12vr2 + vr^) (100) 
3 

The polygamma part is a bit more complicated as it has poles not only at 
the kinematical and dynamical locations but additionally alo q = i{Q + 2n), 



^^In contrast to the Konishi case we have not symmetrized the polygamma part of scalar 
integral. This choice leads to a slightly different configuration of the poles. 



36 



n> 0. The contribution of its residues are 



Res Y^'\q) 

q=i{Q+2n) 



^ n\-l + Av?){n + QY{-1 + 4(n + Qf)) 



n=l 



This can be summed over n resulting in polygamma functions. We combine 
this result with the contribution of the remaining poles obtained by taking 



residues of Y^^ (q) 



+ sP{q,Q) 



at the kinematical and dynamical 



poles. The final expression can be summed over Q by the methods of Ap- 
pendix B. The result is simply 



256 



_7,^(_12 + 7r^)^ + 54(4-3C(3))C(3)) 



(101) 



When we combine this result with the result of the rational part we can 
observe that the even C part cancels out and we arrive at the subleading 
wrapping correction 



Ai^o) = -128(12C(3)2 + 32C(3) + 40C(5) - 105C(7)) 



(102) 



We checked that the dynamical residue contributions of the full integrand 
cancel when summed over Q as in the Konishi case providing another con- 
vincing support for our considerations. 



B Summation of terms containing polygamma 
functions 

Here we present a method which enables one to calculate sums containing 
polygamma functions. A typical sum emerging during the 5-loop Konishi 
computation is of the form: 

oo 

5](m) ^ J2 R{Q) i^^'^\Q) m > 0, (103) 

Q=l 
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where R{x) is a rational function of x and ip^'^\x) is the mth polygamma 
function given by the definition 

^(-^\x) = ^-^^, m>l, (104) 

with ip(x) = ip^^'^^x) being the digamma function ip(x) = ^ J^^^ • The 
evaluation of sum fllOSp goes as follows: R{Q) is decomposed as a sum of two 
terms R{Q) = Ro{Q) + Ri{Q) , where Ro{Q) contains the sum of pure inverse 
power terms of the partial fraction decomposition of R{Q) (i.e. Ro{Q) = 
+ + --O' while Ri{Q) contains the rest. In this case the sum fllOSp is 
decomposed into 2 parts as well: 

= Sj,'") + S!™\ (105) 

where 

oo 

Si") = 5^i?,(Q)^(-)(g) m>0, a = 0,1. (106) 

Q=l 

Using a series representation for the polygamma functions the sum Sq™''' can 
be evaluated directly applying the method sketched in section 6. So, hereafter 
we concentrate on the calculation of the sums They are evaluated 

by transforming them into integral expressions calculable with the help of 
Mathematica. The transformation is as follows: the polygamma functions 
are represented by their appropriate integral representations 

oo 

/j^m p—xt 
dt m>l, (107) 

1 - e * 



oo 

^{x) ^ V^W(x) = j dt{^^- ■ (108) 



and the function Ri{x) is represented as the Laplace transform of its inverse 
Laplace transform 

oo 

R^{x) = j dte-""' CR^\t), (109) 
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where CRi^(t) stands for the inverse Laplace transform of Ri{x) given by 
the formula 

r]+ioo 

CR-\t) = — [ dse^'RAs), (110) 
2m J 

rj—ioo 

with T) being an arbitrary positive constant chosen so that the contour of 
integration lies to the right of all singularities in Ri{s). Due to the structural 
difference between integral representations (11071) and (11081) one has to make a 
distinction between the cases m > 1 and m = and consider them separately, 
m > 1 case: 

Using integral representations (I107P and (I109p . S^™"* takes the form 



oo 



^(rn) = ^ / dte-'^'CR-{\t) (-1)"^+! / dt 

Q=i 



1 - e-*' 

which after evaluating the simple geometric sum in Q becomes 

oo oo 



)(e*+*' - 1) 



Then using the identity 



oo 







, , t'"" , , Li^+i(e-*) - C(l + H , , 



where Li„(a;) is the rath polylogarithm function, the integral with respect to 
i! can be evaluated and finally a single integral remains 

oo 

e!™) = (-1)™ / dt CR^\t) r(m + 1) "5^^ + ""^ . (112) 



During the Konishi computation m took the values of 1, 2 and 3 and in all 
cases emerging during the calculations, the integrals (I112p could be evaluated 
by Mathematica. 
m = case 
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Representing Ri (x) as the Laplace transform of its inverse Laplace transform, 
and using integral representation fllOSp for ip^^^x) the sum takes the form: 



oo oo 

Sf) = ^ J dte-'^'CR^\t) J dt' 



e *' e 



f 1 - e-*' 

Again evaluating the simple geometric sum in Q, this can be recast as 

oo oo 

e-* 1 



J dtCR^\t) j dt' i^— 



t' (e* - 1) (1 - e-*') (e*+*' - 1] 
Now exploiting the integral formula 

^dt' 1 ^ = lE + \og{l-e~^) 

\t'{e^-l) (1 - e-*') (e*+*' - 1); 1 - e* ' ^ ' 



with being the Euler's constant, the integral with respect to t' can be 
evaluated and one ends up with an expression containing a single integral 

oo 

Sf) = / dtCR-At) + (114) 
J 1 - e* 



During Konishi computations this formula made it possible to evaluate sums 
containing il)'^^\x) functions. 

Finally, we just note that in the case of m = the sum Eq°^ can also be 
evaluated by the application of formula flll4p . 
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